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Axion like particles (ALPs) are quite generic in many scenarios for physics beyond the Standard 
Model, they are pseudoscalar Nambu-Goldstone bosons, and appear once any global U{1) symmetry 
is broken spontaneously. The ALPs can gain mass from various non-perturbative quantum effects, 
such as anomalies or instantons. ALPs can couple to the matter sector incluidng a scalar condensate 
such as inflaton or moduli field via derivative interactions, which are suppressed by the axion decay 
constant, f^ • Although weakly interacting, the ALPs can be produced abundantly from the coherent 
oscillations of a homogeneous condensate. In this paper we will study such a scenario where the 
ALPs can be produced abundantly, and in some cases can even overdose the Universe via odd 
and even dimensional operators, as long as fx/^i ^ 1? where 4>i denotes the initial amplitude of 
the coherent oscillations of the scalar condensate, cj). We will briefly mention how such dangerous 
overproduction would affect dark matter and dark radiation abundances in the Universe. 


I. INTRODUCTION 

Inflation [T] is extremely successful paradigm for cre¬ 
ating seed perturbations for the large scale structure in 
the Universe [2]. One of the post-inflationary challenges 
is to recreate the Universe we live in - filled with the 
Standard Model degrees of freedom and the dark matter 
particles [3]. The energy extraction from inflaton energy 
density to matter degrees of freedom is non-trivial, and it 
all depends on the couplings, i.e how inflaton couples to 
other matter degrees of freedom, for a review, see SIS]. 

Similar to the inflaton, there could be many light mod¬ 
uli fields present in the early Universe, whose potentials 
are absolutely flat. The moduli fields are natural conse¬ 
quence of any low energy string theory, which appear in 
4 dimensions with almost flat potential [6]. In any super- 
symmetric theory there would be flat directions, which 
are built upon F and D flatness condition, see HI- Typ- 
ically, the moduli could be made heavy by stringy non- 
perturbative effects [8], but still there could be one light 
moduli with a mass around the scale of low supersymme¬ 
try (SUSY) breaking. If the low scale SUSY breaking is 
at O(TeV), then the moduli would typically take a mass 
of that order. Light moduli could be displaced during 
inflation by virtue of quantum fluctuations of order the 
Hubble expansion rate during inflation, see [H HO]. In 
which case, the moduli lighter than the Hubble expan¬ 
sion rate during inflation could be settled with a large 
vacuum expectation value (VEV), i.e. of the order of 
Mp ^ 2.4 X 10^^ GeV. The moduli behaves like a con¬ 
densate [ 7 ] , and begins its coherent oscillations when the 
Hubble expansion rate of the Universe drops to the mass 
of the moduli Q 


^ In principle one should check the flatness condition of the mod¬ 
uli or flat direction potential by taking various SUSY breaking 
inputs in a time dependent background mi¬ 


lt has been recognised, sometime ago, that the energy 
transfer from the inflaton/moduli field to its decay prod¬ 
ucts need not be perturbative. The perturbative decay 
is typically slow, and the time scale is determined by the 
total decay width of the inflaton/moduli. On the other 
hand, there is an interesting possibility of draining the 
energy density of inflaton/moduli via non-perturbative 
effects in a time dependent background - known as pre¬ 
heating m, and [II El. 

The particle production during preheating is indeed 
non-adiabatic, the effective time dependent frequency 
does not evolve adiabatically. The violation of adia- 
baticity condition is the key for particle creation M- 
Another important non-perturbative production mecha¬ 
nism, dubbed as tachyonic preheating, takes place when 
the time dependent frequency becomes imaginary which 
happens for instance when symmetry is broken sponta¬ 
neously [laiiE], or when the trilinear coupling comes to 
dominate the interaction with the inflaton/moduli [17]. 
The process of particle creation, soon becomes non-linear 
due to the mode-mode interactions, which makes the 
problem soon intractable. The epoch of preheating, or 
explosive particle production stops when the back reac¬ 
tion is taken into account [I11IIH1II9]. Nevertheless, dur¬ 
ing the first few oscillations, as long as the energy density 
in the decayed quanta is still below the energy density 
stored in the inflaton/moduli sector, we can trust the 
analytical treatment, see for instance [20] . 

During preheating one can excite scalars [HHni, 
fermions [5T], higher spins [25] and gauge fields [21- 
1^ . Preheating can also generate isocurvature perturba¬ 
tions 12EI, and non-Gaussian signatures |2ll|2H|. One can 
also excite fairly long-lived dark matter particles [HEg. 
In fact, the large abundance of dark matter particle cre¬ 
ated at the time of preheating is the real threat for cos¬ 
mology. 

In this paper we are keen to understand the cosmo¬ 
logical consequences of inflaton or moduli interaction via 
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derivative couplings. For simplicity we will not distin¬ 
guish between inflaton or moduli field, we will call the 
held to be a generic sealar eondensate without loss of 
generality, which oseillates coherently when the scalar 
condensate mass becomes equal to that of the Hubble 
expansion rate. 

One such example of derivative interaction is demon¬ 
strated by axion-like particles (ALPs) dehned here as y, 
for a review, see m- The ALP appears due to the spon¬ 
taneous breaking of a global U{1) symmetry. At low 
energies below the symmetry breaking scale, a Nambu- 
Goldstone boson arises as massless excitations of the an¬ 
gular part of the complex scalar held, which is an absolute 
gauge singlet. The ALPs can couple to the scalar held 
via dimensional 5 and higher dimensional operators. The 
operators are suppressed by the axion decay constant, 
fx ^ 4 dimensions. Typically could be large 

close to the GUT scale or the string scale, or even the 
Planck scale. We will keep as a free parameter. 

Typically, the axion such as in the case of QGD receives 
its mass from instanton j3Th33] , in the string theory con¬ 
text the anomalous U{1) gets broken non-perturbatively 
by the world-sheet instantons, and its mass could be large 
or small depending on the action of the instanton |34] . 

ALPs can act as a perfect dark matter candidate, due 
to its weak interaction with baryonic matter held j35h37] . 
On the other hand, axion produced from the decay of a 
scalar condensate can give rise to a non-thermal distri¬ 
bution of dark matter and dark radiation, depending on 
the mass of ALPs. Of course, one has to make sure that 
ALPs are produced perturbatively/non-perturbatively to 
match the current observations from the Big Bang Nu¬ 
cleosynthesis (BBN) [38], and the current Planck data on 
dark matter and dark radiation abundance [2| 

The goal of this paper is to demonstrate that: 


The paper is structured as follows - in section [IT] we 
will discuss the ALP’s interaction with a homo gene ous 
scalar condensate, i.e. inhaton/moduli. In section [nH we 
study the ALP scalar interaction via dimensional 5 op¬ 
erator, and study preheating of ALPs without and with 
expansion ehects. We will also discuss two regimes: when 
^ 0{Mp) and when < Mp. In section IV, we will 
discuss the ALP’s interaction via dimensional 6 operator 
in two regimes: when ~ 0{Mp) and when < Mp. 
First we will discuss the production of ALPs, without 
expansion, for light, and heavy, > m^, 

ALPs. We will repeat the computation by taking into 
account of the expansion effect. Finally, we will conclude 
in section |V| by alluding to some of the dark matter and 
dark radiation bounds, whose detailed scans will be pro¬ 
vided in a separate publication. 


II. INTERACTION AND EQUATIONS OF 
MOTION 

Let us consider the following Lagrangian with deriva¬ 
tive couplings of the form|^ 

^ 

+ -^(jy^d^xd^X + • • • , ( 1 ) 

where 0 is a scalar condensate, which sets its coherent 
oscillations around the minimum of a quadratic potential 
governed by its mass 

Let us denote y as an ALP with a mass m^, which 
couples to the 0 held with derivative interactions, and 
is the axion decay constant. Here we will consider 
two scenarios: 


the ALPs ean he produeed abundantly from the eoher- 
ent oseillations of a sealar eondensate for a wide range 
of masses, i.e. heavy or light ALPs, as long as f^ < ^i, 
where denotes the initial amplitude of the eoherent 
oseillations of a sealar eondensate, i.e. inflaton/moduli. 


• when the axion decay constant is ^ ^ Mp. 

• when the axion decay constant is /^ < ^ Mp, 

where we have assumed that the condensate is oscillating 
with an initial amplitude ^ 0(Mp)[^ Assuming that 


For the purpose of illustration, we will hx the initial 
value of the amplitude of the oscillations of the scalar 
condensate to be Mp, without loss of generality. What 
matters here is the ratio: fx/^i- The rate at which the 
particle production will take place would be sensitive to 
this value, as we shall demonstrate in this paper. The 
results of this paper will have an important consequence 
for string theory motivated ALPs 013911111112]. 


^ String theory motivated models of inflation face a stern challenge 
from the axion over production from numerous hidden sectors, 
first raised in this paper m- The question then arises - how 
to arrange various mass scales in string compactiflcation models, 
such that a successful reheating without axion domination can 
take place eq]. 


^ For the purpose of illustration, we will consider only dimen¬ 
sional 5 and 6 operators, in fact one can discuss any higher order 
operators, whose behaviour will be similar to the one discussed 
below. 

One can relax this condition and have oscillations around any 
kind of potential, but for many examples, such as moduli os¬ 
cillations, it is a fairly good approximation to consider coherent 
oscillations around a quadratic potential. A small deviation from 
a quadratic potential would not affect the physical outcome dras¬ 
tically, see m- There are exceptions, when the scalar conden¬ 
sate oscillates around a displaced minimum, in which case there 
would be a kinematical blocking due to finite VEV effect, which 
might hamper preheating la. and reheating both M- 
^ Any light condensate, lighter than the Hubble expansion rate 
during inflation can obtain large VEVs during inflation and act as 
a condensate. In principle, we will be able to relax this condition 
also. As we shall see what matters is the ratio: where 

is the initial amplitude of coherent oscillations. 
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(j){t) is a homogeneous condensate, its equation of motion 
in a Friedmann-Robertson-Walker (FRW) background is 
determined by: 

0 + 3H(j) + = 0 , (2) 

where ’dot’ denotes derivative w.r.t. physical time, t. A 
typical condensate can in principle attain a large ampli¬ 
tude during its coherent oscillations. The solution of the 
above equation is given by 


0 (t) ^(t) sin(m 0 t) 

(3) 

^ Al « 

m^t mci)t 

(4) 


operator and higher order operators are negligible. Both, 
a(t), and the Cc;|(t), can be then approximated by 


— 3H T 


Hfx 


~ 3iJ 


1 + V/x 

cos{m^t)/f^ 
1 + $(t) sm{m^t)/f^ 


(9) 


and 




1 + (t)/fy^ 

a? 1 + i>(i) sm{m^t) /’ 


( 10 ) 


Note that during the coherent oscillations, the average 
equation of state is that of matter dominated epoch, the 
Hubble expansion rate is given by: a{t)/a{t) = H(t) ^ 
(2/3t). When the Hubble parameter, H(t) ^ rricfy^ the 
condensate begins coherent oscillations. The equation of 
motion for y field is then given by : 


d^d^x- 


3H 


{d°<p/fx + 4>d°4>/f^y 

{l + ^/f^ + p/2m 


dox 


+ 


mz 


(l + 0//^ + 0V2/2) 


X = 0. 


(5) 


Expanding y in terms of the momentum modes, the 


above Eq. ([^ can be rewritten as: 


Xk + a{t)xk + ‘^l{t)xk = 0, 

( 6 ) 

where 


m ozJU^ , Hfx + (t>Hfl 
a(t =3iL(t + , ^2/or2’ 

1 + P/Jx + 0 /2/x 

(7) 

and 


a 2 + l + <^/P + <(.7/r 

( 8 ) 


are the damping coefficient, and the effective frequency 
squared, respectively. Throughout this paper we will 
follow the positive energy initial solution, i.e Xk — 
i^2ujk. Let us now study Eq. (§) for various sce¬ 
narios. 


III. ANALYSIS OF A DIMENSIONAL 5 
OPERATOR 

In this section, we will consider dimensional 5 opera¬ 
tor, i.e. (j){dx)‘^/fx^ will assume dimensional 6 


respectively. We now proceed to discuss the evolution of 
the dynamical system for two cases: ~ Mp, and when 

and < Mp. 


A. Axion decay constant of order Planck scale: 

- 0{Mp) 


When ^i/fx < where 
approximated by: 


^ Mp^ Eq. © can be 


(y,{t^ — 3H -|- 


fx 


1 - 


3H + rriff) cos(m0t) ~ sin(m0t) 


( 11 ) 


Similarly Eq. (10) can be approximated by: 


^ I 1 - ^ 


P 2/1 

- < 1 - 


fx, 

. 

sir 

Jx 


{m^t) 


( 12 ) 


Since fx ^ 1? damping coefficient can be approxi¬ 
mated by: a [2 Mp cos(m 0 t)//^]/t, during coherent 
oscillations of 0, which is always positive. Eor the same 
reason, uj‘1 is always positive. Hence, one would naively 
think that there cannot be tachyonic instabilities, which 
are triggered by a negative squared frequency [mis]. 

We will see later that the presence of an oscillat¬ 
ing damping coefficient leads to tachyonic instabilities 
However, as we shall see now, these tachyonic instabil- 
ties do not lead to significant enhancement in Xk when 
> <I>i. Nevertheless, the evolution of the time de¬ 
pendent frequency, cj/c, may become non-adiabatic dur¬ 
ing parts of (p oscillations, leading to narrow parametric 
resonances m- We will investigate this possibility in 
coming section. 


An alternative description with canonically normalised x field is 
derived in an Appendix B. 


^ This can be easil y see n by means of redefinition of Xk: see for 
instance Eqs.( |40j )-( |43| > for the dimensional 6 operator, a similar 
transformation holds for the dimensional 5 operator. 
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B. Excitation of ALPs when ~ 4>i ~ 0(Mp) 


Let us first ignore the damping term and the expan¬ 
sion rate of the Universe by setting a = 0 and a = 1. 
Introducing the following dimensionless measure of time: 
rrifpt = 2z ^ 'K one can recast Eq. into the form of 
Mathieu equation m- 


Xk + [^ - 2qcos{2z)]xk = 0 , 


with 


A = 


4(/c^ + m^) 




(13) 


(14) 


and 


Q = 


4> ml 

9_^ 

fxml 


(15) 


where ’prime’ denotes derivative w.r.t. 2 :. The quanta 
of the X field are produced when the evolution of the 
the time dependent frequency becomes non-adiabatic m 

[T4] . i.e. 


2 

k 


uj- 




2 /x 


k‘^ 




3/2 


>1, (16) 


where \(p\ < Now we will consider two scenarios: 


Case -1: When k < m^: 

For modes with k < the adiabaticity condition, 
see Eq. (16), is violated when 


Hence, one can estimate the maximum enhance¬ 
ment in the energy density of a par ticular momen¬ 
tum mode, pk (X \xk\‘^ (see Eq. |A3[ ) p| during each 
oscillation of (j) by 6 pk ^ exp{4pmAt) ^ exp{^/3)^ 
which means that the energy density of the modes 
defined by the mass and momentum range above 
increases at most by a factor of 6 per each oscilla¬ 
tion of <p. 


• Case -2: When k > m^: 

For modes with k > m^, the adiabaticity condition 
is violated, i.e. uJk/^k — I 5 when 



1 ^ ml 


1/3 


(18) 


which implies that — < Therefore, modes 

with momenta ^ < 7 ^ can be excited. 

Here again, q < |(^)^ < and the allowed pa¬ 
rameter space is in the range: 2 {f^/^)q ^ A < 
( 2 g)^/^. Hence, A > 2 g, and the system never 
hits the tachyonic instability of the Mathieu equa¬ 
tion [isms]. 

Again for 4>i < < 1.24>i, the dynamical system 

hits the first resonance band twice per each oscil¬ 
lation of 0 for a time period At < each. The 

characteristic exponent is also < \/3/4. Hence, 
one would expect that the energy density of the 
modes defined by the mass and momentum range 
above to increase at most by a factor of 6 per each 
oscillation of (j). 


!!k < 14. 

~ 2 /x ’ 


(17) 


which implies that q < |(^)^ < in which case 
the dynamical system effectively enters the narrow 
parametric resonance regime, if A>i-q [nun]. 
Eq. ( |T7| also defines the relation between A and q, 

Az^XfxMQ- 

For > <l>i, the system never enters any tachy¬ 
onic instability nani], since A > 2 g in this case. 
However for 4>i < < l-24>i, it hits the first reso¬ 

nance band leading to the parametric excitation of 
ALPs with masses in the range [4 + < 

m^ < ^m(f) and momenta k < m^. This occurs 
twice per each oscillation of (j) for a time interval 
At < each. 

During these time periods, Xk grows as Xk ~ 
exp(/im 0 t), where p being the characteristic Flo- 
quet exponent [45], which for the first resonance 
band is given hy p = A^ — [q — 1)^. For the 
case at hand, p < \/3/4. 


• By taking into account of expansion: 

Let us now include the effect of expansion. In this 
case the amplitude of both 0 and Xk oscillations die 
out with time. In fact, 4> drops to roughly 0.1 of its 
initial value after the fist oscillation. Consequently, 
q drops by a factor of 1000 and A drops by at least 
a factor of 100 . As a result, the system is pushed 
out of the instability band, and the real part of 
characteristic exponent becomes 0 , even before the 
end of the first oscillation of <p. This render the non- 
perturbative production of ALPs in this scenario 
utterly inefficient. 

To conclude this section, one can safely conclude that 
for > 4>i, there is no significant amplification of Xk- 
In other words, there is no non-perturbative particle pro¬ 
duction of X quanta in this case. The decay from 0 con¬ 
densate to X quanta will be primarily perturbative, and 


® See the definitions of Pfc(t), and the energy density of x quanta, 
Px{t): ill appendix [a] 
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it is given by: 


xx) = 


1 rnl 

64tt /2 


(19) 


which is typical of any slow decay of an inflaton/moduli 
via higher dimensional operator. 


C. Axion decay constant smaller than Planck scale: 

Mp: Ghost like ALPs draining energy 


/^=Mp/25, ni^= 3m(f„ k«;m^ 
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FIG. 1. The above figure shows the time evolution of otiz) jmcf, 
(black curve) and (red curve) for = Mpj^'b. 

The evolution of a{z)lm(f) becomes negative roughly half the 
time period of os cilla tions, and becomes singular at the times 
defined by Eqs. ( [^ and this results in tachyonic in¬ 

stabilities. For k < the evolution of uj\{z)lvr?^ becomes 
negative for less than half the time period of oscillations as 
long as ^//x > 1, see the upper panel, whereas for k > 
the evolution uj‘l{z)lm‘^ becomes negative only in the vicinity 
of singular points, see the lower panel. Here time is measured 
in units the corresponding number of (j) oscillations is 

shown on the upper horizontal axis. 

Now, let us consider the case where ^i//x ^ 1, where 


^ Mp. For convenience, let us rewrite Eq. in 
terms of the dimensionless measure of time, 2 ) = 


where 


and 


Xk 


a{z) 

TTlff) 




2 Xk b , 
mi 


rr(^) — 3H T 


TO 0 ^COs{z)/P 
l + ^sm{z)/P ’ 



mt 


1 + ^sm{z)/fx 


( 20 ) 


( 21 ) 


( 22 ) 


Both a(z) and u)l(z) can be negative leading to an ex- 
ponential growth of and therefore copious particle 
production, see Fig. (IH. In fact, both a{z) and 
become singular, when (j){t) = —fx^ 

4"^ = sin“i(/x/$j) -h (2r - l)7r 

-/,,/i>(^*,,) + (2r-l)7r, (23) 

and 

= -sin-i(/,,/$,) + 2r7r 

-/x/^(^*,r) + 2r7r, (24) 

where r = 1,2 ,... denotes the r*** oscillation of ^(t), and 




(4,-1)-, (25) 


refers to the mid of the r^^ time interval ,where is 
tachyonic, see Fig. Q 

The number of singular points a{z) and oo‘^{z) pass 
through depends on the ratio of /x/^i, i.e. smaller 
this ratio is, the more singular points a{z) and ujI{z) 

(t) fr) 

will have. Between the singular points z^ ^ and Zj , 
0 (t) < —/x, which results in a ghost-like behaviour of 
the X field, as the kinetic term becomes negative between 
these singular points. This can be easily seen from the y 
effective Lagrangian, which in this case is given by: 


^x-\0- + 4>{t)/fx)i9xf - 


In fact, the singular points arise when the kinetic term 
changes form positive to negative, and vice versa. At 
the singular points the coefficient of the kinetic term be¬ 
comes zero. As a consequence, (Qy)^, and hence y^ be¬ 
come infinite. Moreover, when 0 < —/x, < 0 become 

tachyonic, see Eq. (22), leading to an exponential pro¬ 
duction of ALPs with mass ^ and momenta 
k < m^. This can be understood from a very rough solu¬ 
tion, Xk cx: exp(|u;/c|/m 0 Az), during the tachyonic phase. 


® Here and through out this paper, we fixed the initial time, i.e. 
the time at which the scalar condensate starts oscillating, to 7r/2. 
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Also, the parametric excitation of ALP is possible as the 
evolution of in time become non-adiabatic, i.e. 


Wjk 

$ TO0 

COs(2() 


fx 

(1 + ^sm{z)/f^y/^ 


>1, (26) 


for k < m^/\l + ^ sm{z)/^ and 




^ m: 




|cos( 2 ;)| 


P 2 k (1 + ^sin(^)//^)2 


> 1 . 


(27) 


for k > m^/|l + ^ sin(^)//^|^/^. 

Nevertheless, one would anticipate that energy of the 
0 condensate can be drained even before the system ac¬ 
tually hits the first singular point signalling the rise of 
the ghost-like behaviour of the y field. To see this, we 
will examine the behaviour of Xk in the vicinity of the 
first singular point, located at tt + /^/^. Expand¬ 
ing a and to first order around that point, Eq. (20) 
becomes: 


earlier. This renders the problem non-linear even before 
the first singular point, whose analysis goes beyond the 
scope of the current paper. 


f^=Mpl25, 



^//x 






X/C = 0, 
(28) 


where for simplicity the effect of expansion was ignored. 
Eq. (28) admits the following solution: 


Xk^e 


— ikz 






+C2 U[l,l,-2i — ^ l + 


/x, 




])} 


(29) 


where ^ | — f ’ Ln{x) is the Laguerre polynomial, 

and U (a, 6, x) is the confluent hypergeometric function of 
the second kind (also known as the Krummer’s function 
of the second kind) [46]. 

Erom Eq. (29), one can see that Xk is divergent at 


the first singular point, because U{a,b^x) diverges log¬ 
arithmically as its argument, x, approaches zero i.e. as 
z ^ IT fx/^- Hence, ALPs with arbitrary masses and 
momenta can be excited in the vicinity of the first singu¬ 
lar point. 

In particular, for ALP with a given mass, momenta in 
the range: 0 < k ^ (4>//^)(m^/m0) are favoured for y 
production. However, by getting closer to the first sin¬ 
gular point, higher momentum modes can be excited; in 
principal modes with infinite momenta can be excited at 
the exact singular point. However, one does have to get 
much closer to the singular point to drain the energy from 
the condensate. See Eig. ([^, where we have shown that 
the energy density stored in the excitations of a light y 
field, px^ is becoming comparable to that of the conden¬ 
sate even before the first singular point. The definition 


of the energy density of the y field is defined by Eq. (A5), 
see appendix. 

Eor heavy ALPs, with ^ m^, the y energy den¬ 
sity becomes comparable to the scalar condensate even 


FIG. 2. The evolution of Pxl{p<t),iCL ^) during part of the 1st 
oscillation of cj) field. 


The above picture will be the case for any odd¬ 
dimensional operator, such as dimensions 5, 7, 9, • • •. 
Eor even dimensional operators, such as 6, 8, 10, • • •, 
the ghost problem will not arise and we will be able to 
discuss the ALPs dynamics safely. 

We now conclude our study of dimensional 5 operator 
by briefly mentioning that it is indeed possible to excite 
copious ALPs with a wide range of masses and momenta 
during the first oscillation of the cj) condensate. Eollow- 
ing the evolution through out will be a daunting task 
analytically, since the energy density of the ALP very 
soon becomes comparable to that of the energy density 
of the condensate, besides the existence of formidable sin¬ 
gularities. Eurthermore, ghost-like behaviour of the ALP 
arises due to the kinetic term oscillating to negative val¬ 
ues. This does not happen for ALP interaction with even 
dimensional operators, which we will study in the next 
section. 


IV. COPIOUS PRODUCTION OF ALPS VIA 
DIMENSION 6 OPERATOR 


Let us assume that the ALP is coupled to the conden¬ 
sate only via dimension 6 operator, then the damping 
coefficient, a{t), can be approximated by: 


OL{t^ — 3H T 




4>4>/fx 

1 + <^ 2 / 2/2 

sin(2m^t)/2/2 
1 + i>2 sm^{m^t)/2p ’ 


( 30 ) 
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similarly, can be approximated by: 




/ / 

1 + 4>‘^/Vl 

K 

1 + ^2 sin^(m0t)/2/2 


(31) 


Figs. ^ and Q show plots of a/rrifp (black curve) and 
^k/'^c!) (blue curve for k and red curve for k 

m^) for ^ and = ^i/25, respectively. We 

now proceed to discuss the evolution of the dynamical 
system for two cases: and /^ < 


A. Axion decay constant of order Planck scale: 

- 0{Mp) 



z=m(f,t 


When ^i/fx ^ where ^ Mp^ the damping coef¬ 
ficient can be approximated by: 




$2 / $2 \ 

~ 3J7 + rritf,-^ sm{2m^t) I 1 - ^ sin (m^t) 1 . 

^Jx \ ^Jx / 

Similarly, Co’^(t) can be approximated by: 


UJu ^ 


n2 


■ 


1 - 




(32) 


+ ^ sin2(m,^t)^ . (33) 


For ^/fx ^ 1, Of — 2t ^ + ^j( 2 m 0 t^/^) ^ sin( 2 m 0 t) is 
always positive, which can be also seen from Fig. . For 
the same reason, uj‘1 is always positive. However, due to 
the presence of an oscillating damping coefficient, tachy- 
onic instabilities appear for short time intervals while (j) 
is oscillating. This can be easily seen from Eqs.([4Q|)-([4^. 
However, as will see below, these tachyonic instabilties do 
not lead to significant enhancement in Xk when 
Moreover as 0 oscillates, the evolution of ujk in time may 
become non-adiabatic giving rise to parametric instabil¬ 
ities which may lead to exponential enhancement in Xk- 
We will investigate this possibility below. 

In order to investigate the existence of such parametric 
resonance, let us for the time being ignore the damping 
term and the momentum redshift, i.e. by setting a = 0 
and a = 1. Introducing the dimensionless measure of 
time, z = 7 r /2 + m 0 t, Eq. can be put into the form of 
Mathieu equation m 


where 


Xk + [^ - 2qcos{2z)]xk = 0 , 

(34) 

k‘^ mA 

/ $2 \ 

(35) 

A = —^ -|-^ 1 




(36) 


FIG. 3. The above figures show the time evolution of ajmcp 
(black curve), and cukl'n^cf) for k ^ (blue curve) and 
k ^ (red curve) for ~ l.lMp. For k, uJk does 

not oscillate and redshifts with the expansion of the Universe, 
whereas for k <C uJk oscillates with double the frequency 
of 0 field, and moves toward its maximum value, rux, as the 
amplitude of the oscillations dies out with cosmic expansion. 
Here time is measured in units the corresponding num¬ 

ber of oscillations is shown on the upper horizontal axis. 


Again, the quanta of the field Xk are produced when the 
adiabaticity condition gets violated, i.e. 


coi 


$2 


TO^m^l sm{2m^t)\ 


4/2 


+ m^{l - sin^{m^t)} 


ZI2 


> 1 . 


(37) 


Let us consider two scenarios: 


• For modes with k < m^: the adiabaticity condi¬ 
tion is violated, when 


~ 4/2 ■ 


(38) 


Hence, q = 1, which implies that the 

dynamical system enters the resonance regime, if 
A > 1 — q. However, here A < < I — q. 

Hence, the dynamical system never hits the reso¬ 
nance bands of the Mathieu equation. Moreover, 
since A > 2 g, the system does not enter the tachy¬ 
onic instability region. Therefore, there is no signif¬ 
icant non-perturbative particle production in this 
case. 


• For modes with k > nix' adiabaticity condi¬ 

tion is violated for 


< 






-I 1/3 


( 39 ) 



















which implies that ^ As a result, q = 

^Jx 

< 1 and ^ < 1-q. Again, 

A is never below 2 q. Therefore, the dynamical sys¬ 
tem never hits the instability bands of the Mathieu 
equation. 

Including the expansion effect and the damping term, 
which is always positive in this case, will make things 
even worse, as far as non-adiabatic excitation of ALPs is 
concerned. Therefore, one can safely conclude that for 
fx ^ fhere is no significant amplification of Xk- In 
this case, the production of the ALPs will be primarily a 
perturbative process. 


B. Axion decay constant smaller than Planck scale: 

^ ~ Mp 

Now, let us consider the case where ^i/fx ^ with 
^ Mp. In this case, the damping coefficient, a, oscil¬ 
lates with a large amplitude until ^ drops roughly below 
/^, see Fig. 0. This is actually crucial for particle cre¬ 
ation. As we shall see, this gives rise to tachyonic insta¬ 
bility, and parametric resonance, leading to exponential 
particle production. The details will be discussed in the 
following subsections. For convenience, let us introduce 


fx=Mpl25 



FIG. 4. Same as Fig. but for = Mp/25. Note that 
the damping coefficient, a, now oscillates to negative values 
twice per each oscillation of (j) field. 

the following transformation: 


Xk = A'/cC 




M'/fx 

1+02/2/2 


dz 


As a result, Eq. (20) becomes 


A' 


+ —A' + ^ Afc=0, 


(40) 


(41) 


with 


^ p rn'^- ^( 00 " + 0 '^) - §^rnd/''^' 

1 + /V2/I 


+ - 


4/^ 


/^2 


(1 + (Z72/|)^’ 

being square of the frequency of Ak oscillation. 


( 42 ) 


1. Non-perturbative excitations of ALPs, ignoring 
cosmic expansion 


We first start by discussing the non-perturbative exci¬ 
tation of X quanta while ignoring the effect of expansion. 
In this case, Cbl becomes 


Col = k^ + 


cos(22;) 

1 + ^ sin^(2:) 


(1 + ^ sin^iz))^ 


(43) 


Moreover, the exponential part of Eq. (40) can now be 
easily evaluated to hold 


exp 


1 c ; 

2 / l + 4^/2f} 



<^2 

27J 


^1 + ^[1 -cos(2+] 


(44) 


This term clearly exhibits an oscillatory behaviour with 
period tt, measured in units of Hence, it does 

not lead to any enhancement of x/c- Moreover upon the 
inclusion of the expansion effect, ^ dies out with time 
and hence such oscillatory behaviour goes away in a few 
oscillation of 0. It is worth noting that the above term 
is independent of both and k. We now proceed to 
discuss non-perturbative production mechanisms for two 
mass regimes: (1) light ALPs with mass ^ and 
(2) heavy ALPs with > m^. 


a. Narrow resonance production for light 
ALPs, with rrix 

We can see from Eq. ([4^ , see also Eig. ([^, that for 
modes with k < (^k becomes tachyonic dur¬ 

ing short time intervals, around z = nir, n = 0,l,2,---. 
More precisely, Cj^ becomes negative during the time in¬ 
tervals, Az = zp — Z-^ where 


z± nn A 


k $ 


$2 


1/2 


for k > m^, whereas for k m^, we have 


z± ~ n7r± 
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This leads to tachyonic excitation of y modes with mo¬ 
menta in the range 0 < k < The up¬ 

per bound on k can be easily seen from Eq. ( [4^ , be¬ 
cause (jj‘1 is positive for higher momenta. On the other 


f^=Mpl25, k=3m0, m^«:k 



FIG. 5. The time evolution of for k — and ^ 

k, where the effect of expansion is ignored. The Cj‘1 becomes 
negative for two short periods of time during each oscillation 
of 4>. Here time is measured in units the corresponding 

number of cj) oscillations is shown on the upper horizontal axis. 

hand, from Fig. ([^, one can see that Cj‘1 sharply changes 
around 2 ; = riTT, before it becomes negative. Hence, 
one would expect violation of the adiabaticity condition 
around the time intervals where uj^ < 0. Indeed, this is 
the case provided that k < More percisely, 

the adiabaticity gets violated during the time intervals, 
Az = , where 


nTT ± 


^4^ 

fc 3 $2 


1/5 


for k > m. 0 , while for k <C m^, 

5±:^n7r±[4/2/$2]V''- 

This leads to the parametric excitation of y quanta with 
momenta, k < 

Recall that becomes tachyonic during the time in¬ 
tervals = 2 :+ — around 2 : = n 7 r, where 2 ;± n 7 r ± 

_ ^]i /2 foj. f. > and 2;± ~ mr 

for k rriff). Hence, the parametric enhancement of Xk 
takes place only for very short tim^eriods around the 
tachyonic dimples in see Fig. In other words, 
both production mechanisms go side by side. 

In fact, both the instabilities discussed above are 
caused by the oscillating damping coefficient, Eq. ( [3Q| ), 
that oscillates with a large amplitude, see Eig. @ , which 
gives rise to deep tachyonic dimples in Moreover, 
parametric instabilities occurs on the sides of those tachy¬ 
onic dimples before Cu‘^ becomes negative. These insta¬ 
bilities lead to exponential production of ALPs. 


It is worth noting here that even when ff>i ^ fx^ 
these tachyonic dimples will be present when k^m^ < 
m 0 ff>/v/ 2 Mp, albeit being way shallower. Hence, one 
would expect that all modes with /c, m^, can be ex¬ 

cited. However, the production of such modes is insignifi¬ 
cant as the enhancement in Xk after each tachyonic phase, 
which is roughly given by 5xk ^ 

at most a factor of 2 per each oscillation of 0 ; conse¬ 
quently, pk gets enhanced by at most a factor of 4 per 
each oscillation of 0. Eor the dimensional 5 operator 
discussed earlier, p^ gets enhanced by roughly the same 
factor per each oscillation of (j). However, this enhance¬ 
ment will go away upon the inclusion of the expansion 
effect. 


f^=Mpl25, k=3m0, 

123456789 10 



FIG. 6. The above figure shows time evolution of pk for 
k = 3m(f) and m^ <C k, where the effect of expansion is ne¬ 
glected. Here time is measured in units m^^; the correspond¬ 
ing number of 0 oscillations is shown on the upper horizontal 
axis. 


In Fig. (|^, we ha ve shown the time evolution of 
defined by Eq. (A3), for = Mp/25, ^ and 

k = 3 m 0 , while ignoring the effect of expansion. One 
can see that pk exhibits an oscillatory behaviour, except 
for short time intervals, where the aforementioned insta¬ 
bilities resides. During these time intervals, pk increases 
rapidly to stabilize at a larger value after each period of 
instability. 

Eor the parameters given above, pk grows roughly by 
a factor of 2000 per each oscillation of 0 , which is a very 
high production rate. 


b. Non-perturbative production of heavy 
ALPs, with 

The perturbative production of ALPs with > 
7710/2 is kinematically forbidden. However, the produc¬ 
tion of such heavy particles can still take place dur¬ 
ing the coherent oscillations of (p condensate via non- 
perturbative effects. 

Similar to the case of light ALPs, becomes tachyonic 
for short time periods, around z = nir, 77 = 0,1,2,---, 
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provided that < vrid)^/\f2fy. This occurs twice 

per each oscillation of 0, see Eq. ( [^ and Fig. 0 . More 
precisely, < 0 during the time intervals, = 2 :+ — 
around ^ = n 7 r, where 


= riTT ± 




1/2 


for k < m^, and 


z± = niT ^ 


'^0 v^/x 

~Y ^ 


^2 


1/2 


for k > m^. During these time periods, tackyoinc exci¬ 
tation of X quanta with momenta in the range 0 < /c < 
^/\/ 2 /^ takes place. 


f^=Mpl25, m^=3m0, k«;m^ 



FIG. 7. Same as Fig. 0 but for = 3m(/> and k ^ 
without taking into account of expansion. 


In Fig. Q, we have shown the time evolution of Cj‘1 
for = Mpl2b^ = SrUffy and k rUff) while ignoring 
the effect of expansion. Similar to the case < /c, one 
can see that for > k, uj^ changes sharply around the 
time intervals of tachyonic instabilities. Again one would 
expect violation of the adiabaticity condition during the 
time intervals where changes sharply. Indeed this is 
the case; the evolution of uj^ becomes non-adiabatic for 
short time intervals, around the tachyonic dimples resid¬ 
ing in the vicinity of z = n 7 r, see Fig. <§■ This leads 
to the parametric excitation of heavy y quanta with mo¬ 
menta k < m^. Moreover, analogous to the case of light 
ALPs, tachyonic and parametric production of particles 
with < k < m 0 ^/\/ 2 /^ is also possible. 

Therefore, during the short time intervals of paramet¬ 
ric and tachyonic instabilities, heavy ALPs with mo¬ 
menta in the range 0 < k ^ m 0 ^/\/ 2 /^ can be excited. 
Again, very similar to the case of light ALPs, one can eas¬ 
ily see from Eq. (30) and Fig. @ that the instabilities 
discussed above are induced by the oscillating damping 
coefficient. 





FIG. 8 . The evolution of \ujk/^k\ i^i time for = 3m(f) and 
k ^ ?Ti 0 , where the effect of expansion is ignored. One can see 
easily that the adiabaticity gets violated around the intervals 
where becomes tachyonic. This takes place twice per each 
oscillation of cj). Here the time is measured in units the 

corresponding number of 0 oscillations is shown on the upper 
horizontal axis. 


f^=Mpl25, m^=3mfp, 



FIG. 9. Same as Fig. but for — 3m^ and k ^ 
without taking into account of expansion. 


In Fig(|^, we show the time evolution of pk for = 
Sniff), k and = ^i/25 while ignoring the effect 

of expansion. One can easily see that pk gets largely 
enhanced at each interval of instability which happens 
twice per each oscillation of (j). However comparing this 
to Fig([^, one can see clearly that the production of ALPs 
with momenta, k > m^,m 0 is way more efficient than 
those with momenta, k < m^. 

Let ns now identify which of the non-perturbative pro¬ 
duction mechanisms discussed above is dominant. For 
this purpose, we consider the two cases: 

• When < k In this case the dimples in 

uj‘1, see Fig. are almost completely on the neg¬ 
ative side, in which case the tachyonic production 
is expected to be significant. Moreover, the para- 
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f^=Mpl25, k=10 m^<k 
123456789 10 



FIG. 10. The evolution of pk in time for — Mpl2^ and 
< k rricf) where the effect of expansion is included. Here 
the time is measured in units the corresponding number 

of 0 oscillations is shown on the upper horizontal axis. 


metric excitation in this case is expected to be in¬ 
significant as the violation of adiabaticity condition 
is expected to be moderate due to the smallness of 
the part of the dimple in residing on the positive 
side. In other words, for < k the non- 

pertnrbative production of y particles mostly take 
place via tachyonic resonance. Fig. (10) shows the 
evolution of pk in time for < k where 

the effect of expansion is ignored. 


• When < k: in this case parts of those 

dimples, on whose sides the adiabaticity condition 
is violated, reside on the positive side allowing a 
room for parametric particle production. One can 
easily see from Fig. (| 6 ) that the non-perturbative 
production of y particles in this case is way more 
efficient. 


dimples in to become less and less deeper with time, 
see Figs. ( pT] ) and ( p^ . Therefore, one would expect 
that the non-perturbative excitation of y quanta will 
seize to be significant after a few oscillations of (j). 

a. Non-perturbative production of Light 
ALPs, with 

For light ALPs with we saw earlier that the 

non-perturbative production of ALPs is dominated by 
modes with momenta in the range ^ k < m 0 ff>/\/ 2 /^. 
Thus, except for the dimples in a;|, it behaves as k"^ jo? ^ 
for which case it redshifts to lower values with the dim¬ 
ples keep being present, albeit becoming shallower, see 
Fig. ([n]); as a result the violation of adiabaticity due 
to the sharp change in on the sides of those dimples 
becomes more and more moderate with time. Hence de¬ 
pending on the scale /^, the non-perturbative excitations 
of y quanta will become less and less efficient with time, 
until it is completely taken over by the expansion of the 
Universe. Nevertheless, as we shall see, we would end up 
exciting a large abundance of ALPs even in this case. 


/^=Mp/25, k=3m0, m^^k 



Therefore, one may conclude that non-perturbative 
production of light ALPs is dominant in the momentum 
range ^ k ^ m 0 ff>/\/ 2 /^ with parametric resonance 
being the dominant production mechanism. This is also 
the case for Heavy ALPs, > rricf), which can be un¬ 
derstood by comparing the pk for k > m^. Fig. (|^, to 
that for k < m^. Fig. §• Clearly, pk in the later case is 
way less than for the former. 


2. Non-perturbative excitations of ALPs, with 
expansion 

Let us now include the effect of expansion. This will 
result in lowering the production rate in two ways. First, 
the energy density of ALPs, pk^ with a given momenta 
dilutes with the expansion, which becomes obvious at 
late times where the excitation of Xk quanta is less 
significant. Second, the amplitude of (j) oscillations dies 
out with expansion, oc causing the tachyonic 


FIG. 11. Same as Fig. §, but now the cosmic expansion is 
included. One can see that the tachyonic instabilities become 
dominant within few oscillations of field. The corresponding 
number of 4> oscillations is shown on the upper horizontal axis. 


In Fig. we have shown the evolution of pk for 

^ rriff) and k = by taking into 

account the effect of expansion. One can see that the non- 
perturbative production of ALPs becomes insignificant 
roughly at the end of the third oscillation of (j). 

Next, we calculat e the energy density of the y parti¬ 
cles, defined by Eq. (|A5|). For = ff>i/25 and k, 

we have shown the time evolution of p^j{p(f)ja~^) in 
Fig. (13). In this case roughly 10% of the condensate’s 
energy get gets transferred to the excitations of the y 
field. However, when < ff>i/27, most of the energy 
stored in the (j> field can be drained in which case the 
problem becomes nonlinear. 
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f^=Mpl25, m^<^k 

123456789 10 



FIG. 12. The above figure shows the evolution of pk in time, 
where the cosmic expansion is included. Here time is mea¬ 
sured in units the corresponding number of cj) oscilla¬ 

tions is shown on the upper horizontal axis. 

/^=Mp/25, mfj) 



FIG. 13. The above hgure shows the evolution of Px/P(j)j 
in time for rux ^ where rric^ = 10^^ GeV. Here time is 

measured in units the corresponding number of <f> oscil¬ 

lations is shown on the upper horizontal axis. This analysis 
includes cosmic expansion. 


f^=Mpl25, k<^m^ 



FIG. 14. Same as Fig. 0- but now cosmic expansion effect is 
included. One can see that the tachyonic instabilities become 
important within few oscillations of cj). The corresponding 
number of cj) oscillations is shown on the upper horizontal 
axis. 



FIG. 15. Same as Fig. but now the cosmic expansion 
is taken into account. The evolution of ujk seizes to be non- 
adiabatic after a few oscillations of cj). The corresponding 
number of cj) oscillations is shown on the upper horizontal 
axis. 


h. Non-perturbative production of heavy 
ALPs, with rrix ^ 

let us now consider the effect of expansion on the non- 
perturbative production of heavy ALPs. We saw that the 
non-pertrbative excitation of modes with k, < rrix is 
way less efficient than for modes with < rrix ^ 

Figs. and (§. Including the expansion effect, the dif¬ 
ference in efficiency between the two cases becomes even 
more pronounced for the following reason. As the ampli¬ 
tude of (j) oscillation decays with time, will have dif¬ 
ferent behaviour for k > rUx and k < rrix- For the former 
case, uj‘^ tends toward smaller values with the dimples in 
uj‘1 keep being present, albeit becoming shallower and the 
violation of adiabaticity condition on the sides of those 
dimples becomes milder with time, see Fig 0- 


On the other hand, for k < m^, tends toward larger 
positive values, i.e toward its maximum value which 
is m^, with the dimples becoming shallower with time 
and even disappearing in a couple oscillations of 0, see 
Fig ( p^ . This behaviour can be also seen from Fig. (^. 
With the disappearance of the dimples in d;^, the adia¬ 
baticity condition is no more violated which can be seen 
from Fig. (15). Therefore, one would expect that the 
non-perturbative excitation of modes with k > rrix 
longer that that for modes with k < rrix- 

This render the non-perturbative particle production 
in the later case less efficient which can be also seen from 
Figs. (12) and (16) that show the time evolution of pk 
for k > rUx and k < m^, respectively, where = ^i/25. 
For the later case, the particle production seizes to be 
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f^=Mpl25, m^=3m(i), 

123456789 10 



FIG. 16. The above figure shows the time evolution of pk- 
As expected, the excitation of such heavy ALPs takes place 
during the first few oscillations of cj). Here time is measured 
in units the corresponding number of (j) oscillations is 

shown on the upper horizontal axis. This analysis includes 
cosmic expansion. 


efficient roughly after the first (p oscillation, whereas it 
keeps being significant roughly till the end of the third p 
oscillation for the former case. 


f^=Mpl25, 

123456789 10 



FIG. 17. The above figure shows the time evolution of px/pcj),! 
for rrix = where = 10^^ GeV. Here time is measured 
in units the corresponding number of p oscillations is 

shown on the upper horizontal axis. This analysis includes 
cosmic expansion. 


Now we proceed to calculate the p^ for tak¬ 

ing into account the effect of expansion. In Fig. (17), we 
have shown the relative energy density of ALP to that of 
p particles, respectively. It is evident that particle cre¬ 
ation is way less efficient for heavy ALPs as compared to 
light ones. This is because of the smallness of the win¬ 
dow for the more efficient non-perturbative production of 
modes with momenta k > (i.e < k < 
with rrix ^ in contrast to the case < rUff). 

Therefore, in order for the non-perturbative produc¬ 
tion of heavy ALPs to be efficient, a smaller fx/^i is re¬ 



FIG. 18. The above figure shows Px/p(t>,i as a function of 
ratio fx/^i after 5 p oscillations for = 10^° GeV (dashed 
curves) and = 10^^ GeV (solid curves). The black (red) 
curve corresponds to light (heavy) ALPs, respectively. This 
analysis includes cosmic expansion. 


quired. For instance, when > ^i/47 and 
most the condensate’s energy density can be transferred 
to the X field, within few oscillations of the scalar con¬ 
densate, p. Heavier ALPs requires even smaller <C 4>i. 
In Fig. (18), we have shown a plot of the energy density 
of X field to the redshifted initial energy density of the 
scalar condensate as a function of the ratio fx/^i after 
5 p oscillations. The black (red) curve corresponds to 
light (heavy) ALPs, respectively. One can easily see the 
difference in efficiency between the two cases. It is also 
clear that the ratio fx/^i being ^ 1 is crucial for efficient 
particle creation including cosmic expansion. 


V. CONCLUSION 

In this paper, we have considered the non-perturbative 
production of ALPs during the coherent oscillations of 
a scalar condensate with mass m^. The ALPs couple 
only derivatively to other fields including inffaton or any 
moduli field. We focused on the couplings, dimensional 
5: p{dxY/fx^ dimensional 6 operators: 0^(9x)^//^- 

In both cases, when > 4>i, there is no signifi¬ 
cant non-perturbative production of ALPs. Neverthe¬ 
less, when fx 4>i, the non-perturbative excitation of 
ALPs is significant. However, in case of dimensional 5 
operator, p{dx)‘^/fx^ singularities in both the damping 
coefficient, a and the frequency squared, arise when 
p = —fx- These singular points occur when the kinetic 
term changes from being positive to negative and vice 
versa giving rise to an unphysical ghost-like behaviour of 
X between these singular points, i.e. when 0 < — /^. 

The energy of the scalar condensate can still be drained 
even before the dynamical system hits the first singular 
point rendering the problem non-linear before the first 
singular point. This will be the case if any coupling with 
odd powers of p dominates the interaction with ALPs. 

On the contrary, such ghost-like behaviour does not 
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arise in even dimension interactions, In this 

case, ALPs with masses in the range 0 — m04>/A/2/^ 
can be copiously produced via non-perturbative effects, 
provided that <C 4>i. If is sufficiently below, 4>i 
(/^ < 4>i/27 for ^ and even smaller values for 
heavy ALPs, > m^, for instance < 4>i/47 for 
= Sniffs), most of the energy density of the scalar 
condensate would be transferred to the y excitations. Ir¬ 
respective of inffaton or moduli as a scalar condensate, 
this will definitely cause an alarming problem for ALP- 
dominated Universe, as it would lead to significant de¬ 
parture from the Standard Model cosmology and could 
spoil the success of BBN, unless the ALP is unstable and 
decay before the era of BBN. These constraints will be 
discussed in a separate publication. 

Moreover, if y is stable enough to be the dark mat¬ 
ter, one can show that ALPs with a wide mass range can 
overdose the universe. Only in a very fine-tuned region 
of the parameter space can one match the observed DM 
abundance. Furthermore, the recent Planck constraint 
on the number of relativistic species, A'eff = 3.15lo.4o E]? 
which places an upper bound on the abundance of dark 
radiation, < 2.92 x 10“^, can be used to fur¬ 

ther constrain the parameter space toward the lower mass 
end mm?]. 

In principle, ALPs could also couple to other fields. 
For instance, it could couple to the Higgs via the opera- 
tor; or to a gauge boson via the operator: 

I which in the ALP’s rest frame leads to a 
decay rate: F ^ ALPs could also couple to 

fermionic fields via the operator: fx' 

decay rate of ALPs will be then given by: F ^ m^m ^//^, 
in the y rest frame. The perturbative decay rates above, 
being already slow, are even slower for lighter ALPs, 
rUff), which are relativistic, due to the smallness 
of their mass and the time dilation effect when the decay 
rates are transformed from the particles rest frame to the 
cosmic comoving frame. 

This can be used to further constrain the ALP’s pa¬ 
rameter space. For instance, if ALP constitute a com¬ 
ponent or the whole of DM, their mean life time should 
be longer than the age of the universe which places an 
upper bound of 40 MeV (/^/lO^^GeV)^/^ on the mass 
for ALPs. Heavier ALPs would decay at earlier stages 
and hence, their mean life time and consequently, their 
mass would be subject to several constraints such as from 
CMB, BBN, and other astrophysical observations. 

In view of our results, it would be indeed interest¬ 
ing to revisit some of the earlier computations on ax- 
ion production in the string theory context, where most 
of the analysis was done in a perturbative decay of the 
infiaton/moduli into axions [39fi42] . Our current study 
highlights that by taking into account of various non- 
perturbative effects, the Universe after string theory mo¬ 
tivated infiationary models could be filled abundantly 
with non-thermal, both heavy and light axions. Al¬ 
though, we have discussed the ALP production for a sin¬ 
gle scalar condensate, with multi-scalar condensate, the 


situation may get even worst, such as in the case of as¬ 
sisted inflation [48]. 

In future, it would be also interesting to introduce, sys¬ 
tematically, higher order interactions and sum them up 
in a natural way such that the theory of axion interaction 
becomes an infinite derivative theory [49| , very similar to 
the gravitational interaction [50] . 

Before we conclude, the most important aspect of our 
analysis suggests that the ALPs can be overproduced in 
the early Universe for ^ 4>i during the coherent oscil¬ 
lations of the infiaton or a moduli condensate, the sum¬ 
mary plot is Fig.[^ which shows how the energy density 
in the ALPs for a wide range of masses can overdose the 
Universe for fx/^i < 1- I^i P^st much of the attention 
has been paid to understand perturbative production of 
ALPs, but the non-perturbative production of ALPs is 
perhaps the most dangerous channel to drain the infiaton 
or moduli energy density. In light of this discussion, it 
is important that we reconsider some of the analysis of 
axion production in the early Universe. 
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Appendix A: ALP’s energy density 

The stress-energy tensor for the Lagrangian under con¬ 
sideration, Eq. 0, is given by 

- (^1 + ^ + daXdpX 

(Al) 


Hence, the energy density of the y field is then given by 



(A2) 


where the prefactor (l + 0//^ + 0^/2/^) exhibits an os¬ 
cillating behaviour with amplitude decaying with the ex¬ 
pansion of the Universe. When the amplitude of (j) oscil¬ 
lation, 4>, drops below /^, this prefactor can be ignored. 
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Furthermore, the energy density for a given momentum 
mode would then be given by 


Pk (l+y^ + 2/2 


1 


+ 7;<eff\Xk\^ 


(A3) 


where ujk.effi denoting the effective frequency of Xk os¬ 
cillation, is defined by 


1/2 


Wfe,e// = [^ + 


1+-A+ 


(A4) 


fx 2/2 


For simplicity of the notation, we will call ujk.eff^ just 
ujk‘ In terms of the ALP’s energy density is given by 

Px = ^ J dk k'^Pk ■ (A5) 


Appendix B: Canonically normalised Lagrangian 


Instead of our approach, one could have first canon¬ 
ically normalized the kinetic term in the Lagrangian 
Eq. 0 . which could be easily achieved by following field 
redefinition: 


X ^ (1 + 0//x + / /v) 


X- 


The Lagrangian then reads 


^2f (t> I 0^ W/2 fdn(f> I 

1^„_ , 1 2i2 ^ 1 — + + 


2 2,1 A ^ 1 ,1 (^ + ^) ddi^4>. 2 


C= -d^ct>d^ct>--m^cP +-d^xd^^X-i^rn^X +/ ^ ^ 


xd^X+a' 




2 2-x' '2 (1 + A + ^) + m + A + u h 

(Bl) 


r 


Assuming that 0 remains a homogeneous condensate, one 
can then derive the following equation of motion for the 
X field: 


d^dfj,x + ^Hdox + 


/dg0 0do0 w d 0 ^ 

t /x fx fx 

4 n ^ 

+ 2/2 J 


+ 


- 

2 V 


(f)dod^(f) 

P 


!) \ S/f ( do(f) I (/)do(/) \ 

P 2 V ; 


(1 + ^ + 4) 


X = 0. 


(B2) 


r 


Expanding x in terms of the momentum modes, Eq. (B2) 
can be rewritten as 


Xk + 3iJ (t)xfe + Wfe {t)xk = 0, 


(B3) 


where 


— “2 + 


7 9 1 { (fp _i_ 4>4> \ (p \ 3H ( (p (p(p\ 

e “ 2iA + - —Ia + TFi 




1 ^A + Ipfx + /p 

”^4 n + A. + Ai')2 

+ A ^ 2/d 


(B4) 


When the interaction with the scalar condensate is dom¬ 


inated by the dimensional 5 coupling, Eq. (B4) can be 


approximated by 


‘^k{k) — -y + 


_0_/ 0 \2 

2f„ 2/’,W \2f„! 


(1+^) (1+^)2 


(B5) 


On the other hand, when the dimension 6 operator domi¬ 
nates the interaction with the scalar condensate, Eq. (B4) 
can be approximated by 


^l{t) ^ 




ml - \{^ + ^)- If# 


P 


2/5 


( 00 ^2 
I 2 f 2 ) 


(1 + 


(i + &^ 


(B6) 

Eq. (B3) with frequency uJk defined by Eq. (B6) for 


the dimensional 6 operator becomes exactly the same as 
Eq. (41) with ujk defined by Eq. (42), along with ^ 

One can also show the same happens for dimensional 5 
and all higher order operators. Both the approaches dis¬ 
cussed in this paper, i.e. with and without canonical 
kinetic term, would yield similar numerical result for the 
energy density, p^, stored in the ALPs. This has been 
numerically verified. 
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